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ABSTRACT 

We study the global dynamics of advection-dominated accretion flows (ADAFs) with 
magnetically driven outflows. A fraction of gases in the accretion flow is accelerated into the 
outflows, which leads to decreasing of the mass accretion rate in the accretion flow towards 
the black hole. We find that the r-dependent mass accretion rate is close to a power-law one, 
mcxr % as assumed in the advection-dominated inflow-outflow solution (ADIOS), in the outer 
region of the ADAF, while it deviates significantly from the power-law r-dependent accretion 
rate in the inner region of the ADAF. It is found that the structure of the ADAF is significantly 
changed in the presence of the outflows. The temperatures of the ions and electrons in the 
ADAF decreases in the presence of outflows, as a fraction of gravitational power released in 
the ADAF is tapped to accelerate the outflows. 
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1 INTRODUCTION 

It is widely believed that many astrophysical objects are powered 
by mass accretion on to black holes. The standard geometrically 
thin, optically thick accretion disc model can successfully explain 
most of the observ ational features in active g alactic nuclei (AGN) 
, and X-ray binaries jshakura & Sunvaevlll973h . In the standard thin 
■ model, the motion of the matter in the accretion disc is nearly Ke- 
I plerian, and the gravitational energy released in the disc is radi- 
1 ated away locally. An alternative accretion disc model, namely, the 
advection-dominated accretion flow (ADAF) mode l, was suggested 
for the black holes accreting at very low rates ( Ilchimanj|l977l ; 
iNaravan & "Yil ll994l) . In the ADAF model, only a small fraction 
of the gravitational energy released in the accretion flow is radiated 
away due to inefficient cooling, and most of the energy is stored 
in the accretion flow and advected to the black hole. The ADAFs 
are optically thin and hot (comparable with the virial temperature 
of the ga ses in the flows), wh ich radiate mostly in X-ray wave- 
band (see lNaravan & McClintoc k 2008, for a review and references 
therein). This model can successfully explain the main observa- 
tional features of bla ck hole X-ray binaries and l ow-luminosity 
AGN (LLAGN) (e.g . iNarayan & Yilll994 Il995ij ; iGammie et al.l 



llgumenshchev et al."2003': 'McKinnevI 120061 : iGammie et alJI 19991 : 
.Ouataert et a l. 199 9: Yuan et al. 2003). 

iBlandford & BegelmanI d 19991) proposed a self-similar 
advection-dominated inflow-outflow solution (ADIOS) for the 
ADAF with winds. In ADIOS model, the mass accretion rate is no 
longer a constant and is assumed to be a power-law dependence 
of radius (m oc , < s < 1), whic h is an important ingredient 
in most of the follow-up works (e.g. , lOuataert & Naravanll999l : 
lYuan et alj l2003l : IXue & Wanj |20()5|) . Moti vated by the results 
of numerical simulations on accretion discs, Ixie & YuanI j2008h 
investigated the influence of outflows on the accretion flow based 
on a 1.5-dimensional description of the accretion flow. They 
suggested that their solutions can be described by a power-law 
?--dependent mass accretion rate fairly well. 

Magnetic fields are believed in accretion flows, and the 
magnetorotational instabilit y (MRI) provides the s o urce o f vis- 
cosity in accretion flows jsalbus & HawlevI Il99ll . 1 19981) . The 
outflows/jets can be driven by th e large-scale ordered mag netic 
fields threading the accretion disc dBla ndford & Pavri3ll982h . The 
physics of magnetically accelerated outfl ows has been extensi vely 
explo r ed in many previous wo r ks (e.g., Cao & Sprui t 1994: Cag 



1 19991 : lOuataert et all 1 19991 : lYuan et alj l2003l :" lllol |2008|) . As the 
Bernoulli parameter of an ADAF is positive, the ADAF is likely 
to have an outflow, which was confirmed by numerical simula- 
tions and also supported by observations dStone & Pringi3l200ll : 



20021: iKudoh & Shibatj 



Kato, Kudoh. & Shibata 



Naravan. McKinnev. & Farmer 



1995 : iKoide. Shibata. & Kudoh 



2002; lOgilvie & Livic 



20071 : iMcKinnev & Naravad 

2007al lhl). Such outflows/jets magnetically driven from the accre- 
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tion discs provide an ef ficient angula r momentum loss mechanism 
for accretion discs (see ISpniill2008l , for a review and references 
therein). The structure of a standard thin di sc/ ADAF may be altered 
by the magnetically driven outflows (e.g., iLi, Wang. & GMill2008l : 
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In this work, we investigate the global 
structure of an ADAF with magnetically driven outflows/jets. 



2 MODEL 

We consider a steady ADAF with magnetically driven outflows/jets 
surrounding a black hole in this work. 
The continuity equation is 

d 



dR 



(27ri?E«R) + 4nRm^ 



0, 



(1) 



where vr is the radial velocity, E = 2Hp is the surface density of 
the accretion flow, and riiw is the mass loss rate from unit surface 
area of accretion flow. The half-thickness of the disc H is given 
hy H = Cs/^K, and SIk is the Keplerian angular velocity. The 
sound speed Cs = (P/p)^^^, and the total pressure P is the sum 
of the gas pressure and the magnetic pressure: P — Pgas + Pm = 
Pi + Pc + Pm (Pi and Pc are the ion pressure and the electron 
pressure respectively). 

In this work, we adopt the Paczynski-Wiita potential 



GM 



(2) 



to simulate the general relativistic effects of a Schwarzschild black 
hole, where M is the mas s of the black hole, and R ^ = 2GM/(? 
is the gravitational radius l lPaczvhski & Wiitalll98(3h . 
The radial momentum equation is 



1 dp 
I^Ir 



5m = 0, 



(3) 



where Q. is the angular velocity of the accretion flow. The radial 
magnetic force is given by 



5m 



(4) 



27rS ' 

where B^ and B~ are the radial and vertical components of the 
magnetic fields at the disc surface. 

The angular momentum equation reads 

1 d 



0, 



(5) 



dj^R') 1 d a Tm 

where a-viscosity Trip = —aP is adopted l lShakura & SunvaevI 
Il973h . and Tm is the magnetic torque exerted on the accretion flow 
due to the outflows/jets. The outflow is accelerated by the magnetic 
fields threading the rotating accretion disc, and therefore the torque 
Tin can be calculated with 



Tm = 2m^Q.{Ri){R\~ rI), 



(6) 



where rhw is the mass loss rate due to the outflow, i?d is the radius 
of t he footpoin t of the field line, and Ra is the Alfven point (see, 
e.g.. lCaoll20 Q2l. for the details). 

The energy equations for ions and electrons are given by 



2mw£c 
2H 



= 0, 



(7) 



, dsc Pc dp . . + 
P^'^^dR-^dR^-^'' -^'^ + ' 

and 

,dei Pi dp. . ^. + 2rh^ei 

p^'^^dR - - - + + ^ = 0' 

respectively, where the parameter 5 describes the fraction of the 
viscously dissipated energy that goes directly into electrons in the 
accretion flow, and the specific internal energy of electrons and ions 
are given by 



£i = 



7c — 1 pcm-R 

1 kTi 
7i - 1 fumn ' 



(9) 



(10) 



where Tc and Ti are the temperature of electrons and ions respec- 
tively, and the mean molecular weight of the ions and the electrons: 
Pi = 1.23, Pa ~ 1.14 are adopted. The adiabatic indices of the 
electrons and ions, 70 and 71, are given by 



7o 



7i 



1 + 1 



3A'3(l/6'e)+A'l(l/6(e 



4A-2(l/ec) 

3A-3(lM) + Ai(l/ei) 
4A2(1M) 



(11) 



(12) 



where K's are the modified Bessel functions, and the dimension- 
less electron and ion te mperature are defined as: 6c = kTc/ {mcX?) 
and Oi = kTi/{mpC^) jNaravan & ^Il995bl) . The energy dissipa- 
tion rate per unit volume is given by g"*" = —aPRdQ/dR, and 
Qic indicates the energy transfer rate fro m ions to electrons throug h 
Coulomb collisions, which is given by jStepnev & Guilberd[l983f) 



Qic 



3 ruc {kTi - kTc) 

~- 2;;^"^"""^^A'2(l/ec)A2(lM 

2{ec + eif + 1 fec + 



■InA 



(^o + ^i) 



(13) 



where the Coulomb logarithm InA = 20, and is the radiative 
cooling rate consisting of synchrotron, bremsstrahlung, and Comp- 
ton cooling (see, Narayan & Yi 1995b; Manmoto 2000, for details). 

The dynamical properties of magnetically driven outflows/jets 
from an accretion disk can be investigated by solving a set 
of magneto-hydrodynamical (MHD) equations if the magnetic 
field configuration of the disk and suitab le boun dary condi- 
tions at the disk surface are supplied (e.g., ICao & Spruit 1994 



Caol2002l:lKudoh & Shibatall995l:lKoide. Shibata , & Kudoh 199J; 
Kato, Kudoh, & Shibatall2002l ; IOgilvie & Lividll998i.i2001i) . How- 



ever, the generation and maintenanc e of large-scale magnetic fields 
of the disk is still quite unclear (e.g.. lLubow, Papaloizou. & Pringld 
|T994a. b.; Tout & Pringlelll996l) . In this work, we focus on how the 
dynamics of the ADAF is affected by the presence of magneti- 
cally driven outflows/jets. For simplicity, we assume that large- 
scale magnetic field lines thread the accretion disk, and the strength 
of the magnetic fields far from the disc surface along the field line 
to be roughly self-similar: 



Bp{R) ~ B 



pd 



RdJ 



(14) 



where Rd is the radius of the field footpoint at the disc surface, 
Bpd is the strength of the poloidal component of the field at the 
disc surface, Bp{R) is the fie ld strength at R along the field line, 
the self-similar index C ^ 1 dSlandfo rd & Payne" 19 8^, and C = 
4 is ad opted in the calculations of [Lubow, Papaloiz oo, & Pringld 
l ll994bl) . 

The self-similar wind solution derived by 
iLubow. Papaloizou. & Pringld l ll994bl) is only valid for slowly 
moving (non- relativistic) outflows, which was extended for rela- 
tivistic jets by lCaql ll 2Q02| ). In this work, we model the magneti cally 
driven outflows/jets with the approac h adopted in lCad ( l2002h . We 
summarize the model as follows (see lCaoll2002l for the details). 

For a relativistic jet a ccelerated by the magnetic field of the 
disc, the Alfven velocity is l lMichellll969l :lc amenzinal 1 9861) 
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VA = 



(47rpA7j)^''^' 



(15) 



where Bp and pA are the poloidal field strength and the density of 
the outflow/jet at Alfven point, and 7j is the Lorentz factor of the 
bulk motion of the outflows/jets. In this work, all our calculations 
of the outflows/jets are in the special relativistic frame. The Alfven 
velocity va ~ RA^{Rd), where _Ra is the radius of the Alfven 
point along the field line, and Q,{Rd) is the angular velocity of the 
accretion flow at the field footpoint R^. 

The mass and magnetic flux conservation along the field line 
requires 



rriw ^ PAVA 
Bpd 



B^ 



(16) 



where m„ is the mass loss rate in the outflow/jet from unit surface 
area of the disc. 

The final bulk velocity of the magnetically driven outflow/jet 
is ~ Va, so the Lorentz factor of the outflow/jet is 



7j 



(17) 



Combining equations lll4t-l|17|l, the mass loss rate in the out- 
flow/jet from the unit surface area of the disc is 



47rc 



Rd^{Rd) 



(18) 



The origin of the ordered magnetic fields threading 
the disc is still unclear. It was suggested that the mag- 
netic fields can be generated throu g h dy n amo proces ses 



in the disc (e . g.. IShakura & Sunvae' 



evi 



Il99d lArmitagd Tl998f : iRomanova et al 

scale external magnetic fields are transported mward by 
the accretion flow (e.g., Bisnovatvi-Kogan & Ruzmaikid 1 9761 ; 



I973I : iTout & Pringi3 
or the large- 
transported inward by 



iLubow, Papaloizou, & Pringld Il994al : ISpmit & Uzdenskvl l2005h . 
For simplicity, we assume the strength of the large-scale magnetic 
fields threading the disc to be comparable with that of the fields in 
the accretion disc, Bpd — B. The magnetic pressure is convention- 
ally assumed to be proportional to the gas pressure in the accretion 
flow. Thus, we have 



8tt 



(19) 



where /3 is the ratio of the gas pressure to the total pressure, B is 
the strength of the magnetic fields in the accretion flow. 

The magnetically driven outflow is described by the terminal 
velocity of the outflow (Alfven velocity va) when the values of two 
parameters ( and /3 are specified. 



3 RESULTS 

We use the Runge-Kutta method to solve a set of five differential 
equations (I), (3), (5), (7) and (8) for five variables: p, hr, Q, Te and 
Ti with suitable boundary conditions at the outer radius -Rout. In 
our calculations, we adopt the black hole mass M = 10^ Mq for 
a typical AGN. The conventional values of the disc parameters: 
Q = 0.1 and 5 = 0.1, are adopted in all our calculations. The tem- 
perature of the ADAF at th e outer radius is adopte d as described by 
the self-similar solution o f lNaravan&^ ( ll995d) . We use a shoot- 
ing point method in our calculations. Integrating these five equa- 
tions from the outer boundary of the flow at i? = -Rout inwards 



toward the black hole, we can obtain the global structure of the 
accretion flow passing the sonic point smoothly to the black hole 
horizon by tuning the value of radial velocity at -Rout. The outer 
radius of the accretion flow -Rout ~ 5000-Rg is adopted. We find 
that the structure of the ADAF is insensitive to the outer boundary 
conditions. As discussed in Sect.[2] the magnetically driven outflow 
is described by ua, C ™d P- The terminal bulk velocity of the out- 
flow is comparable with the Alfen velocity iia- For any unbounded 
outflows, its bulk velocity should > vk, which implies va ^ vk- 

The mass accretion rate at the outer radius, AI = lO^^'-A/Edd 
(the Eddington rate is defined as MEdd = 1.5 x 10^*i\//MQ 
g s~^), is adopted for the calculations plotted in Figs.[T]-[3] In Fig. 
[T] the global solutions for the ADAFs with outflows are shown with 
different magnetic field strengths and distributions (i.e., different 
values of /3 and Q, in which va = '^k is adopted. In the left panel 
of Fig. [T] the radial velocity and the sound speed as functions of 
radius with different values of /? are plotted for <^ = 1. In the right 
panel of Fig. [T] the mass accretion rates as functions of radius are 
plotted for the global solutions with different values of /3 and i^. 

In Fig.|2l we plot different quantities of ADAF solutions with 
= 1 and 4, respectively, where i^a ~ vk and j3 = 0.9 are 
adopted. As comparison, we also plot the ADAF solutions with- 
out outflows in the same figure. In the calculations, the terminal 
velocity of the outflow is a free parameter. We compare the global 
ADAF solutions with different values of va (i.e., va = vk and 
VA = 2vk) in Fig.[3]with [3 = 0.9 and C = 1. In Figs.lljje] we 
plot the results calculated with a relatively high mass accretion rate, 
mo — 10^^, at the outer radius of the disc. 

The ratio of the bolometric luminosity of accretion disc Lboi 
to kinetic power of outflows Pk.w is calculated with 



-'^bol 



flout (--Rout 

q-'iirRHAR/ I (7j - l)mwC^47ri?d-R.(20) 

-Rin -'-Rin 

We plot the ratio -Lboi/f'k.w as functions of magnetic field strength 
P in Fig.|7] where va ~ Vk is adopted. 



4 DISCUSSION 

A fraction of gases in accretion flow is carried away by the mag- 
netically driven outflow, which leads to mass accretion rate of the 
accretion flow decreasing towards the black hole. In many previ- 
ous works, the outflow is induced by assuming the mass accre- 
tion rate m to be a power-law dependence of radius mocr" (e.g., 
iBlandford & Begelmanlll999l) . In this work, we obtain global so- 
lutions of ADAFs with magnetically driven outflows. Our results 
show that the mass accretion rate rh decreases towards the black 
hole, which is close to a power-law r-dependence at larger radii, 
while it deviates from a power-law r-dependence in the inner re- 
gion of the ADAF close to the black hole (see Figs.[T]and|4ll. The 
large-scale magnetic fields are assumed to thread the accretion flow, 
which are believed to accelerate the outflow from the accretion disc. 
In this work, the magnetic strength at the disc surface is described 
by a parameter /3, which is limited by the gas pressure in the ac- 
cretion flow. For comparison, we also plot power-law r -dependent 
accretion rates in Fig. [T] Our calculations show that the mass loss 
rate in the outflow cannot be very high even if the magnetic fields 
of the ADAF are very strong (see Fig.lTJ, i.e., the mass loss rate 
in the outflow is less than a power-law r-dependent accretion rate 
with s < 0.73. We also calculate the cases with high accretion 



qualitatively similar to those with m = 10 
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For the ADAFs with magnetically driven outflows, the struc- 
ture of the accretion flow is significantly altered in the presence of 
the outflow. We find that both the ion and electron temperatures of 
the ADAF decrease with increasing mass loss rate in the outflow 
(see Figs. |2] and O. This is due to the fact that a fraction of the 
gravitational energy released in the accretion flow is tapped to ac- 
celerate the outflow, which decreases the heating of the ADAF. It is 
easy to understand that the ratio Lhoi/Pk.w increases with decreas- 
ing magnetic field strength (see Fig.|7}. Based on the magnetically 
driven outflow model described in Sect. 2, the kinetic power of the 
outflow Pk,w can be estimated with 



-Pk,w oc -m„vX <x — T— -, (21) 

in the non-relativistic limit. It is found that the kinetic power of the 
outflow decreases with increasing terminal outflow velocity, which 
leads to the structure (e.g., temperatures) of ADAFs is less altered 
by the outflows with higher terminal velocity (see Figs.[3]and[6j. 

In our present calculations, the magnetic field strength is es- 
timated with the gas pressure in the ADAF, which is true for the 
fi elds generated wi th dynamo processes (see, e.g., the discussion 
in lLivio et al In this work, the fields are implicitly assumed 

to be balanced between diffusion and dynamo/inward advection in 



Global dynamics ofADAFs with outflow 5 




10° 10' 10" 10' 10* 10° 10' 10" 10° 10* 



s 




Figure 2. The global structures of ADAFs with magnetically driven outflows. We also plot the structures of ADAFs without outflows for comparison (black 
lines). 



our calculations. The ordered magnetic fields may be maintained 
by the drag-in process in the accretion disc, in which the magnetic 
field strength is determined by the b alance of the advection and dif- 
fusion of the fields in the disc ("e.g.. lLubow. Papaloizou. & Pringld 
Il994ah . The magnetic fields can therefore be stronger than the 
equipartition limit, however, the detailed physics is still quite 
unclear (e.g.. iLubow. Papaloizou^ & Pringlel[l9 94b; Cao & Spruill 
I2OO2I : ISpruit & Uzdenskvl l2005l : iouan & Gammia ,20091) . In this 
case, the value /3 (see Eq.[T9} describing the magnetic field strength 
may vary with radius, which is determined by the balance between 
the diffusion and advection of the fiel ds in the accretion flow (e.g., 
iLubow. Papaloizou. & Pringlelll994bh . The mass loss rate in the 
outflow can be higher than the results presented in this work, if the 
advection of the magnetic fields by the accretion flow is dominant 
over the diffusion in the disc. 
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Figure 5. The same as Fig. [2] but the accretion rate at the outer radius r'rio = 10 ^ , and /3 = 0.925, are adopted. 
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Figure 7. The ratio of the bolometric luminosity of accretion disc to the kinetic power of outflows for different values of (va = "fk is adopted). The black 
lines and red lines are for the accretion rate at the outer radius mo = 10"^ and 10~^ respectively, while the solid lines and dashed lines correspond to f = 1 
and 4 respectively. 



